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Abstract. We find two-sided inequalities for the generalized hypergeometric 
function of the form q+ iF q {— x) with positive parameters restricted by certain 
additional conditions. Both lower and upper bounds agree with the value of 
q+ iF q (— x) at the endpoints of positive semi-axis and are asymptotically precise 
at one of the endpoints. The inequalities are derived from a theorem asserting 
the monotony of the quotient of two generalized hypergeometric functions with 
shifted parameters. The proofs hinge on a generalized Stieltjes representation 
of the generalized hypergeometric function. This representation also provides 
yet another method to deduce the second Thomae relation for 3^2(1) and leads 
to an integral representations of 4^3 (x) in terms of the Appell function F3. In 
the last section of the paper we list some open questions and conjectures. 
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1. Introduction. We will use standard notation for generalized hypergeometric function: 



q+lFq (( a g+l); (bq); %) = q+l F q 



,) 



n=0 



( a l) n(fl2)n ' ' ° (^q-l-l)n n 

(h) n ■ ■ ■ (bg) n n\ 



where (a q ) = a%, 0,2, ■ ■ ■ , a q is the parameter list and (a) n = a{a + 1) • • • (a + n — 1) is the shifted 
factorial. The series converges for \z\ < 1 and can be analytically continued onto the entire complex 
plane cut along [l,oo]. The celebrated continued fraction of Gauss 



2 Fx (a, b; c; -x) 



c — 1 b(c — a)x a(c — b)x 



iFi(a — 1, b; c — 1; — x) c— 1 + 



+ c+1 + + 

(b + n)(c — a + n)x (a + n)(c — b + n)x 



+ c + 2n 

implies on setting a = 1 the classical representation 

c— 1 b(c— l)x (c — b)x 



+ 



c + 2n + 1 



+ 



(2) 



iFx(l,b; c; -x) 



c- 1 + 



+ c + 1 + 



(b + n)(c — 1 + n)x (n + l)(c — b + n)x 



+ 



c + 2n 



+ c + 2n + 1 



+ 



(3) 



institute of Applied Mathematics, Vladivostok, Russia, e-ma,i\:dmkrp@yandex.ru 

2 Voronezh Institute of the Ministry of Internal Affairs of the Russian Federation, e-m&i\:box2008in@gmail. com 



1 



This continued fraction converges and has positive elements when x > and c > b > 1. Hence, its 
even convergents form an increasing sequence and approximate the value of 2Fi(l,b; c; —x) from 
below, while the odd convergents form a decreasing sequence and approximate 2 -Pi (1, b; c; —x) from 
above (see details in [11 ] 1131 [2Tj ) . Taking the first three terms we get: 

1 n /-, ; \ C(C + 1) + (C - b)x 

< 2 F 1 (l,b;c;-x) < / ' \ / < 1. (4) 



l + 6x/c v ' ' ' ' c(c+ 1) + c(6 + l)x 

Gauss derived his continued fraction ([2]) from contiguous relations for the hyper geometric function. 
Another way to explain both the continued fraction (J3j) and the inequality @ is through Euler's 
integral representation 



which shows that 2F1O-, 6; c; —x) is a Stieltjes function, i.e. a function of the form 

1/R 

/(*) = 



d(f>(u) 



1 + xu 







with a bounded, nondecreasing function taking infinitely many values (see [U Chapter 5]). For 
such functions both a continued fraction representation with positive elements and a two-sided 
estimates through Pade approximants are well known (see [4l [5J [TTJ [13]) . Indeed, it is easy to check 
that 

' 1- -x + 0(x 2 ), x^O, 



1 + bx/c c 

c(c + 1) + (c - b)x b ,6(6 + 1) 2 , 3^ 

-1 x + —, -x z + 0(x ,i ), x^O. 



c(c+ 1) + c(b + l)x c c(c + l) 

Hence, the fraction on the left-hand side of (J4j) is the Pade approximant of 2F1O-, 6; c; —x) of order 
[0/1] while the fraction on the right-hand side is the Pade approximation of order [1/1] (see [HE]). 
The sequences of diagonal and paradiagonal one-point Pade approximants form upper and lower 
bounds uniformly converging to the Stieltjes function, i.e. 

[M/M - l] f < [M - 1/M] f < f < [M/M] f < [M + 1/M - l] f (6) 

for each positive integer M. Moreover, these bounds are best possible with respect to the given num- 
ber of power series coefficients [H Theorem 15.2]. Inequality (|3D is a particular case when M = 1. 
These results have been recently generalized to multi-point Pade approximants by Gilewicz, Pindor, 
Telega and Tokarzewski (see [13 EH]). These papers deal, however, with diagonal and superdiagonal 
approximants. The subdiagonal approximants used here are not discussed. The convergents of the 
continued fraction ([2]) were considered by Belevitch in [8]. Explicit expressions for two-point Pade 
approximants for ratios of the Gauss hyper geometric functions, confluent hyper geometric functions 
and Bessel functions are found in [31] . Some applications of Pade approximants to inequalities for 
special functions are discussed in [27] [28] . 

The lower bound in (|4|) is not only asymptotically precise at x = but also agrees with 
2-Fi(l, b; c; —x) at x = 00. One problem with (jH) is that the same, unfortunately, is not true for the 
upper bound, which reduces to the constant (c — 6)/[c(6 + 1)] at x = 00. If we wish to generalize 
PI to q+\F q we are also faced with the problem that neither an analogue of the Gauss continued 
fraction for the general q +iF q is known nor a multiple Euler representation generalizing ([5]) has the 
form of a Stieltjes function. In this paper we solve both problems and find two-sided estimates for 
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f{x) = q+ iF q (l, (a q ); (b q ); —x), where the lower bound is asymptotically precise at x = 0, the upper 
bound is asymptotically precise at x = oo and both agree with the values of f{x) at the endpoints 
of [0, oo]. Next, we extend our inequalities to q+ \F q (o, (a q ); (b q ); —x) for some values of a. These 
results are derived from a somewhat stronger statement about the monotony of a special quotient 
of hypergeometric functions. Our method is based on a generalized Stieltjes representation for 
q+ iF q ((a q+ i); (b q ); —x) but does not utilize the relationship with Pade approximants. In the last 
section of the paper we list some open problems and conjectures. 

Inequalities for general q +iF q are surprisingly rare in the literature. Most important results are 
due to Luke [22], [23] who uses inequalities between diagonal and sub-diagonal Pade approximants 
for (l + z)~P and repeatedly integrates them with respect to beta-distributions or Laguerre distribu- 
tions. See section 2 for some details. Carlson studied in [10] some inequalities for i?-hypergeometric 
function, which can be expressed in terms of a Lauricella Fp function and as such is a generalization 
of 2-Fi to the multivariate case. For 2 Fi(a, b; c; x) his restrictions on the parameters are: c > b > 0, 
x < 1. See section 3 for detailed comparison of our results for 2 F± (a, b; c; x) with those of Carlson. 
Inequalities for q+ iF q (l, (a q ); (b q );x), q > 1 are not considered in Carlson's paper. 

Buschman [9] uses determinant representations to obtain two-sided inequalities for the Gauss 
hypergeometric function 2 Fi (a + n, b; c; x) in terms of 2-P1 (a, b; c; x) for positive parameters and 
x € (0, 1). His results were later improved and extended by Joshi and Arya [THl [19] . 

Inequalities of a different nature for \ p F q \ and "tR.( p F q ) have been obtained by Jahangiri and 
Silvia in [17] for the special case when p F q is subordinate (in the sense of analytic function theory) 
to the linear fractional transformation. 

The papers [J fT2"j [25j [29] (and many more found in the references therein) consider 
inequalities for the Gauss function 2-^1(0, b; c; x), the Kummer function i-F\(a; c; x), the Bessel 
function oF\(a;x) and their ratios for x € (0, 1) (cf. [27]). The emphasis in these papers is on fine 
properties near the singular point x = 1 of 2-P1 (a, b; c; x). Note also that the method for proving the 
monotonicity of ratios based on Lemma 2.1 from [25] cannot be applied to prove our Theorem [TJ 
Since our results are valid also for x < 0, they complement those from [U EJ [TJ [25] for 2-^1- General 
p F q is not discussed in the above papers, except some results in [29] that are applicable when certain 
upper and lower parameters of p F q differ by half-integers. 



2. Generalized Stieltjes representation and inequalities for q +\F q . We begin with a 
representation of ^+1^(2:) as a generalized Stieltjes transform. 



Lemma 1 For Wj^ > ^flfc > 0, k = 1, 2, . . . , q, and \ arg(l + z)\ < n we have 

N 1 



z\=A 



s^- 1 g((a q y,(b q y,s)ds 
{1 + szY 



(7) 



where 



g((a g ); (b q ); s) = J [I - s/(t 2 ■ ■ ■ t^-^ 1 f[ ^"^(l - t^-^ 1 dt 2 --- dt g , 



k=2 



A q {s) = [CI,!]"- 1 n{t 2 ,...,t q : t 2 ---t q > s} 



and 



A n = T 



(bq) 

(a q ),(b q -a q ) 



r(6a)...r(6,) 



r(ai) • • • r(a,)r(6i - ai ) ■ ■ ■ T{b q - a q ) 



(8) 
(9) 
(10) 



The function g({a q ); (b q ); s) is invariant under simultaneous shifts of all parameters: 



g((a q ) + 6; (b q ) + 5; s) = g((a q ); (b q ); s) 



(11) 
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for any complex 5. 

Proof. The multiple Euler integral representation for q+ iF q (z) reads |26} formula 7.2.3(10)]: 



q+\F q 



( K) 






-)- 







(a q ),(b q -a q ) 



[0,1]" 



dh-'-dtg. (12) 



Integration in (fT2j) is over the g-dimensional unit cube, W^k > 3?«fc > for all k = 1, 2, ... , g, and 
| arg(l + z)| < 7r. Formula (fT2j) is obtained by repeated application of the generalized Euler integral 
[21 formula (2.2.2)] 



q+lF q 



(a q ), a q+ i 



T(a q+1 )T(b q - a q+ 



i / 



<zi dt, 



where the last step with q = 1 is the standard Euler integral for 2-^1- 

We make the variable change s = t\t2 ■ • • t q , t\ = s/(t2 ■ ■ ■ t q ), leaving t2, ■ ■ ■ , t q unchanged. The 
Jacobian of such transformation is 



J(s, t2, ■ ■ ■ ,t„ 



t2 ' • • tq 



Hence, we get from (I12|) : 



q+lF q 



a, (a q ) 



z = A Cl 



ds 



(1 + sz) a 



A,(-) 



/ [s/(h ■ ■ ■ t q )r-i[i - s /(t 2 . . . t q )r-^- 1 n t a r\i - t k ) h ~ a ^ 

k=2 t2 "'^ 
[1 - Sf(t 2 ■ ■ ■ tg)}^- 1 f[ t a k ^-\l - t k ) h *- a *- X dt 2 --- dt, 



s° 1-1 ds 
(1 + sz) a 



k=2 



where A g (s) and Aq are defined in ([9]) and (fl~0j) . respectively. Introducing the notation (jHJ) for the 
inner integral we arrive at formula ([7]). The shift invariance (lllh is obvious from the definition of 
9((a Q ); (b q );s). □ 

The representation (JT]) is a key ingredient in the proof of the following theorem. 



Theorem 1 Suppose 5 > 0, > at > 0, A; = 1, . . . , q. The function 
f(a, (a q ); (b q );5;x) = 



q+lFq (([, (Og) + $} (bq) + ~x) 

q+iFq (a, (a q ); (b q ); -x) 



(13) 



is monotone decreasing if a > and monotone increasing if a < /or all x > —1 
Proof. Put 

^ = r 



(bq)+S 

(a q ) + (5, (6, - a q ) 



Then by © and (HU): 



/(cr, (o g ); 



Mo* 



1 „oi+5-l 



(1 + sx)" CT c/((a (? ) + <5; (6,) + 5; s)ds 



Aq Jo s-i-^l + sx)- CT 5 ((a g ); (b q );s)ds 
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_ A S Jo s^-^l + S x)-°g((a q ); (b q ); s)ds 

A o Jo s a ^ 1 (l + sx)-"g((a q );(b q );s)ds ' 

The statement of the theorem is equivalent to saying that f' x (<J-, (a q ); (b q );S;x) < when a > and 
f' x (cr, (ag); (6 ? );<5;x) > when a < 0. Differentiating the definition of / above we see that both 
inequalities are equivalent to the single inequality 

ill l 
f q ( S ) P (s) d s f h ( S ) P (s) d s < J q (s)Hs) P (s)ds J P (s)ds, (14) 



where 

p( s ) = s ^+ 5 - 1 (l + S x)- CT 5 ((ag);(6g 
</<») - A 



1 + 

The function p(s) is positive, while the functions q(s) and h(s) are monotone increasing for fixed 
x > — 1 and < s < 1. Hence, the above inequality is an instance of the Chebyshev inequality [241 
Chapter IX, formula (1.1)]- O 

The value of /(cr, (a g ); (6 g );<5;x) at x = oo can be found using the representation of q+ \F q (z) 
in the neighborhood of the singular point z = oo. Assume for the moment that no numerator 
parameters differ by an integer. Then, according to [26], formula 7.2.3.77], 



lFq [ (M 



(M 



9+1 
fc=l 



«fc) (°?+l) _ «fc 

(bq) - a k 



q+\F q 



l + a k - (b q ),a k 
1 + a k - (a q+ i)' 



1/z) z 



-a k 



(15) 

where the prime at (ctg+i) means that the term a k is excluded from the list. It follows that for 
a < min(ai, fl2j ■ • • , a>q) we will have 



q+\F q 



a, (a q ) 
(b q ) 



(Off) - <*> (M 
(flq),Q> q )-<T 



z~ a + o(z~ a ), z^oo. 



Hence, 



/(cr, (a q );(b q );S;oo) = T 



(a q ), (a q ) + 5 -a, (b q ) - a, (b q ) + 5 
(a q ) - a, (a q ) + S, (b q ), (b q ) + 5 -a 



(16) 



(17) 



When some of the numerator parameters on the left-hand side of (|15p differ by an integer, formula 
(|15[) breaks down and one has to resort to much more involved [261 formula 7.2.3.78]. Formula (|17p. 
however, remains valid by continuity ( p F q /T[(b q )] is an entire function of its parameters - see [261 
7.3.2.8]). 

Formula (|17|) and Theorem [T] imply the following inequalities valid for x > 0: 



/(cr, (o g ); (6g); 5; oo) < /(cr, (o g ); (6,);*; x) < 1 = /(cr, (o g ); 5; 0), 
for positive a and 

1 < /(cr, (o 9 ); (6g);5;x) < f{a, (a q ); (b q );5;oo), 

for negative a. 

Taking 5 = 1 in (|17|) we will have: 

f(a,(a q );(b q );l;o )=f[ b ^^\. 

fJl ai(bi - cr) 



(18) 
(19) 



(20) 



5 



A simple calculation shows that 

<? 

1 - q +iF q (l, (a q ); (b q ); -x) = x q+1 F q (l, (a q ) + 1; (b q ) + 1; -x) JJ(aj/6j). 
It follows that 



i=l 



1 - q+ iF q (l, (a q ); (b q ); -x) _ q+ \F q (l, (a q ) + 1; (b q ) + 1; -x) 

/(l,(a g );(6 g );l;x)J];(oi/^ 



x q+ iF q (l, (a q ); (b q ); -x) 



q+1 F q (l,(a q );(b q );-x) n^A) 



(21) 



i=l 



Combined with (|18|) and (|20|) this gives: 



n 



(«i - 1) l - («?); -«) 



11(6,-1) - * g+1 F g (l,(a 3 ); < II^)- 



< 



Formulas (|18p and (|20p imply that these bounds are best possible and by Theorem [T] the function 
in the middle is monotone. A simple rearrangement of the last formula leads to 



Theorem 2 For bj. > a*. > 1, k = 1, . . . , q, and x > the inequality 

< q+lF q {l, (a q ); (b q ); -x) < 



i+^nLiKA) 

holds true. 

Some comments are here in order. Since, clearly, 



1 + znLiKoi 



(22) 



1 



1 + x Ili=l( a i/bi 



l-xil^- + 0(x 2 ), x^O, 



i=l 



we conclude that the lower bound in (|22p is asymptotically precise at x = and coincides with 
the Pade approximant of order [0/1] to q+ \F q (l, (a q ); (b q ); — x) at zero . Hence, the inequality 
from below in (|22p can be proved by noticing that q+ iF q (l, (a q ); (b q ); —x) is a Stieltjes function by 
Lemma [TJ This allows to relax the restrictions on the parameters to b^ > > 0. Moreover, a 
result of Gilewicz and Magnus [15] (see also [24} Chapter 25, Theorem 2]) implies that the lower 
bound in (1221) is valid for all x > — 1. We further generalize the lower bound in the following 
theorem. 



Theorem 3 For b^ > > 0, k = 1, . . . , q, x > — 1 and a > 1 the inequality 

1 



(1 + xUlM/h)) 

holds true. 

Proof. The proof is an application of Jensen's inequality 



a < q+iFq(a, (a q );(b q );-x) 



(23) 



I b \ 

J f(s)dfi(s x 

a 

b 

fdu 



ftp(f(s))du(t 



< 



b 

fdu 

a 



(24) 
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valid for convex ip, integrable / and non-negative measure [i [24, formula (7.15)]. Take (p(u) = u a , 
<r>l, f(s) = l/(l + sx),dfi(s) = s^- 1 g({a q y,{b q );s)ds. By (7j) 



i 

J f(s)dfj,(s) = ^9+1-^9 (1) (a*); M 





1 1 

j d^ = j s a ^ 1 g((a q );(b q ); S )ds 



^ ^ q+1 F q (l,(a q );(b q );0) = l_ 

A) ^0 
o 



and 



l 

J f(f(s))dfi(s) = ^9+1^ (°"> K); (M; - x ) • 



Hence, <^^ reads 

(g+l-P? (1, K); (69); < g+1^9 (f, (o g ); (6 g ); -x) . 

Combined with the lower bound from (|22p this yields (|23p . The restrictions on the parameters are 
explained in the remark preceding this theorem. □ 

The inequality (|23|) was previously obtained by Luke in [22] using a completely different method 
(see Introduction). His restrictions are bt > dk > 0, k = 1, . . . , q, cr > and a? > so that 
Theorem [3] extends the validity of Luke's inequality to — 1 < x < under the additional restriction 
a > 1. The case < a < 1, —1 < x < remains open. We conjecture that (|23p is still true in 
this case. Another curious lower bound due to Luke [22] valid for bk > Ofe > 0, k = 1, . . . , q, x > 0, 
< cr < 1 is given by 

a+-nL(«^)) < K); (M; 

In a similar fashion but under stronger assumptions we can prove a generalization of the upper 
bound from (1221): 



Theorem 4 For b/. > a/~ > 1, k = 1, . . . , q, x > and < cr < 1 the inequality 

q+1 F q (a, K); (6 9 ); -*) < {1 + xUU[(a '_ _ (25) 

/ioWs trite. 

Proof. Again apply (|24|) but this time with (p(u) = u 1/a , < a < 1, /(s) = 1/(1 + sx) <T and 
cfy/(s) = s ai-1 <7((a g ); (b q );s)ds. This yields 

(g + iF g (cr, (a 9 ); (6 9 ); -x)) 1/a < q +iF q (l, (a q ); (b q ); -x). 

The combination of this inequality with the upper bound from (|22p results in (|25p . □ 
Formula (|16|) and the relation 

1 

x — > 00, 



i+*nu("i -!)/(<>< -1)] 



i5T = ;II^T + W). 



show that the upper bound in (|22p is asymptotically precise at x = 00. This, unfortunately, is not 
true of the upper bound ([25]) . We conjecture an asymptotically precise upper bound for a > 1 in 
the last section of the paper. We also remark that, unlike lower bounds, our upper bounds are very 
different from those of Luke [221 23J. 
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Remark. It is interesting to observe that the proofs of Theorems and H] work for any 
generalized Stieltjes function, i.e. any function of the form 

l/R 

Ux) = 







(i + xty 

with a bounded nondecreasing function (ft taking infinitely many values and a > 0. We have then 

[fi(x)T < [MOT^U^), <r>l, 

and 

fa(x) < [fim^lfl^W, 0<a<l. 

Hence, by combining these inequalities with inequality ([6]) we can "exponentiate" the known in- 
equalities between a Stieltjes function and its Pade approximants. 

3. The case q = 1. In this case we are able to extend the upper bound from (|25p to negative x 
as follows: 

Theorem 5 For b > a + 1 > 1, < a < 1 and — 1 < x < the inequality 

2 F 1 (a, a; b; -x) < - ^ (26) 

(l + 5*1*) 

holds true. 

Proof. The function / defined by (|13|) for q = 1, S = 1 and a = 1 takes the form 

2 Fi(l,a + + l;-x) 



f(l,a;b; l;x) 



2 F 1 (l,a;b; -x) 



The value /(l, a; b; 1; —1) is finite under the assumptions of the theorem and is found by the Gauss 
formula for 2 F\(ai,a 2 ;b; 1): 

fh , , n 2 Fi(l,a + l;6 + l;l) b 
f{l,a;b; 1;-1) = 



2^(1, a; 6; 1) 6-1' 
By Theorem Q] /(l, a; 6; is monotone decreasing for x E (—1,0) and so 

/(l,a;M;0) = l</(l,o;6;l;x) < ^-j = /(l, a; 6; 1; -1). 

According to (|2T|) 

a 1 — 2-^1(1, a; —x) a ^ & 



6 X2-Fi(l, a; fr; — x) 6 \b — 1 
or 

1 1 

— < 2 F 1 (l,a;b; -x) < — — — — . 

1 + ax/b 1 + ax/(b — 1) 

Inequality ([26]) is obtained from this estimate by an application of Jensen's inequality as in the 
proof of Theorem [H □ 

It is interesting to compare these results with those from [10J. Carlson's inequality for the case 
b > a > a > considered here reads 



max 



' l+ '' r ' ' T -,(l+ax/b)^\< 2 F 1 (a,a;b;-x)<(l + xr^ b . (27) 



(1 + x(l - a/b)) b ~ 
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Note first that our condition 6>a>a>0is not more restrictive than b > max(a, a) > 0, since 
one can exchange the roles of a and a. We see that the lower bound in (|27p is an extension (and 
possibly a refinement due to the competitive term under max) of our inequality (|23|) to the values 
of a G (0, 1) for the particular case 9=1. 

To compare the upper bounds we note that for b > a > a > according to ([TBI) 



and clearly 



2 Fi(a,a;b; -x) 



(1 + z) 



x u 



-aa/b 



T(b)T(a-a) 

r(o)r(6 - a) 



+ o(x 



oo , 



1 

P<jafb 



:i+o(i/x)), 



X 



oo. 



Hence, the upper bound in (|27|) is never asymptotically precise, albeit agrees with the value = 
2-Fi(a, b; c; — x) at x = oo. Our bound in ([22]) and the conjectured bound in (|38|) are asymptotically 
precise and so both are better than Carlson's bound at least for large x. The upper bound (|25p 
agrees with the main asymptotic term in order of x but not in the constant, while the upper bound 
in (I27p has the wrong order and so again our bound (|25p is better than Carlson's bound at least 
for large x. 

Finally, for negative x the upper bound in (127h goes to oo as x — > —1, while under our restrictions 
on parameters 2F\(a,a;b; —x) remains bounded and so does our bound in (|25p . Hence our bound 
is guaranteed to be better around x = — 1, while Carlson's bound is more precise when x is close 
to 0. 



3. The case q = 2. In this case we are able to relax the assumptions on the parameters imposed 
in Lemma [T]and Theorems [lj El [3] and [H To this end we give an alternative proof of Lemma [T] which 
also shows that the function g defined by ([8]) is expressed for q = 2 in terms of 2-Fi- Representation 
(|28p below is not new, it is a slightly different form of |26[ formula 2.21.1.26]. However, we include 
a short proof which clarifies the source of restrictions on parameters. 

Lemma 2 Let K(d + e — b — c) > 0, > 0, 3^6 > and \ arg(l — x)\ < it, then 

3 F 2 (a,b,c; d,e;x) 

r(d)T(e) } ^(l - t) d+e - b - c ~ 1 



T{b)T{c)T{d + e-b-c) 



/i (T t) c ~ 
1 J 2 Fi(e -c,d-c;d + e-b-c;l- t)dt. (28) 



Proof. Expand (1 — xt) a on the right-hand side of (|28p into a binomial series and integrate term 
by term applying the generalized Euler integral [21 Theorem 2.2.4] 

l 

U 7-1 (1 _ u y-l-\ Fl {a, (3; 7 ; usjdti = ^^^2^1(0, # u\ x), 
u 

where > > 0, | arg(l — x)\ < tt (for this formula to be valid at x = 1 the additional restriction 
3f?(i> — q — /?) > must be imposed), and the Gauss formula 



F(n ,Hl- r ^- a -^ 
2 F 1 (a,/3,M)- r0y _ a)r> _ /3) , 
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valid under the same restriction. This yields: 
l 







T(d + e-b-c)T(b + k) 
T(d + e- c + k) : 



■Fi (e — c, d — c; d + e — c + k; 1) 



T(d + e-b- c)T(b + k)T(c + k) 
r(e + k)Y(d + k) 



(29) 



which implies (J28J) on substitution into the series. □ 

Remark. Representation (|28p provides another way to derive Thomae's second fundamental re- 
lation for 3 F 2 (a,b,c;d,e;l) [31 formula 3.2(2)]. Indeed, apply the connection formula [2, formula 
(2.3.13)] 

^ / , , , . T(d + e — b — c)T(c — b) , _ . 

2 -Fi(e -c,d-c;d + e- 6-c;l-t)= v — y \ ^ 2 -Pi(d-c,e - c; 6 - c + 1; t) 



+ 



r(d- 6)r(e - 6) 
r(d + e-6-c)r(6-c) +c _ b 



r(d - c)r(e - c) 

on the right-hand side of (f28l) , change t — > 1 — t and apply (f28l) again to get 



* c_ 2Fi(d-b,e-b;c-b+l;t) 



3^2 



a,b,c _ 7rr(d)r(e)r(l - a)r(d + e - b - c) 
d,e J " sin(7r(c - 6))r(6)r(c) 

3F 2 [1 — a, 1 — c, <Z + e — a — 6 — c;d — a — c+l,e — a — c+1] 
r(d - a - c + l)r(e - a - c + l)r(d - 6)r(e - b) 



idem(6; c) 



where the unit argument is omitted for conciseness and idem(6; c) after an expression means that 
the preceding expression is repeated with b and c interchanged. To obtain the form given in [31 
formula 3.2(2)], one needs to apply Thomae's first fundamental relation 



3-^2 



d + e — a — b — c, 1 — a, 1 — c 
d — a — c+1, e — a — c+1 

r(d - a - c + l)r(e - a - c + l)r(fr) 
" r(d + e - a - b - c)T(l + b- a)T(l + b-c) 



iF 2 



b,l + b-e,l + b-d 
1 + b-aA + b-c 



■ (30) 



Using Lemma [2] we can give a version of Theorem Q] for the ratio of 3-F2S under slightly weaker 
restrictions on parameters. 

Theorem 6 Let d + e — b — c>0, c>0, b > 0, 5 > and min(6, c) < min(d, e). Then the 
function 

3 F 2 (a,b + 5,c + 5;d + 5,e + 6;-x) 

/(a, b, c, d, e, d; x) = — — r 

3 F 2 {a,b,c; d, e; -x) 

is monotone decreasing if a > and monotone increasing if a < for all x > —1. 
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Proof. Assume without loss of generality c = min(6, c) (otherwise exchange the roles of b and c). 
This and the hypotheses of the theorem imply that d — c > and e — c > 0. Now follow the proof 
of Theorem [1] to get inequality (I14p with 

S b ~ 1 (l — s \d+e-b-c-l 

p(s) = — 2 F 1 (e - c, d - c; d + e - b - c; 1 - s), 

(1 + xs) a 

J 



s) = s d , h(s) 



1 + xs 

The function p(s) is positive since e > c, d > c and d + e — b — c > 0, while the functions g(s) and 
/i(s) are monotone increasing for fixed x > — 1, < s < 1. The result follows by the Chebyshev 
inequality as before. □ 

Inequality ([22]) takes the form 

' < 3 F 2 (l,b,c;d,e;-x) < — — - — — \-— — —. (31) 



1 + xbc/de ° v ' ' ' ' ' 1 l + a;(6-l)(c-l)/(d-l)(e-l) 

and here is valid under the assumptions of Theorem [6j A particular case of this inequality has been 
used in [20] to obtain an error bound in the asymptotic expansion of the first incomplete elliptic 
integral which leads to very precise two-sided inequalities for this integral. Theorems [3] and H] are 
also valid here under the assumptions of Theorem [6l 

4. The case q = 3. In this section we will show that Lemma [1] leads to new representations of 
4-F3 as a double integral of 2-F1 or as a single integral of the Appell function F3. To this end we need 
to demonstrate that the function $(((23); (63)55) can be expressed in terms of the Appell function 
F3 or as an integral of 2-^1- Indeed, we have by 



4 F 3 (ct, (a 3 ); (63); ~z) 
1 

/„a\ — l f l„ r 
7— / [1 " s/{t 2 t Z )p-^-h a 2 ^- l tf-^-\l - t 2 ) h ^-\l - tsfs-^dhdts, 

(l + sz) a J 

AsO) 

where 

A = r(6!)r(6 2 )r(6 3 ) 

r(oi)r(o2)r(o 3 )r(6i - ai)r(& 2 - a 2 )r(6 3 - o 3 ) ' 

As(s) = {t2,t-z : t 2 t 3 > s,0 < t 2 < 1,0 < t z < 1}. 
The double integral ([8]) here can be written as follows: 

1 1 

g((a 3 );(b 3 );s) = J t^(l - hf^ 1 dt 2 J [t 2 t 3 - s^-l^-b, (1 _ ^a-aa-l^. 

s s/t 2 

Make the change of variables 

X= T ~77^> h = x(l - s/t 2 ) + s/t 2 , 1 -t 3 = (1 - s/t 2 )(l - z), dt 3 = (1 - s/t 2 )dx, 

1 - s/i 2 
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in the inner integral to get 
1 

g((a 3 );(b 3 y,s)= [ (1 - h)^' 1 



dt 2 



x J [t 2 x(l - s/t 2 )] b ^-\x(l - s/t 2 ) + s/t 2 ) a ^({l - s/t 2 ){l - x))^-^- 1 ^ - s/t 2 )da 
o 



s a-3-bi J ^.b 1 +a 2 -ai-a :j -l 



(1 - t 2 ) b2 ~ a2 ~\l - s/t 2 ) bl+bs - ai - a3 - l dt 2 ' ' 



bi— oi— 1 



(1-x) 



-a 3 -X 



(l + x(t 2 /s-l)) b ^ 



■dx 



-ai)r(63-a 3 ) 6l 
r(6i + 6 3 - a\ - 03) 



1 

< y ^+02-«l-«8-l(l _ 4 2 ) 6 2-^-l(- 1 _ g / t2 



,6i+68-ai-09-i_ P . I h-a 3 ,bx-ax 
bx + h-ax- a 3 



2-f 7 ! 



1 - I di 2 



r(6i - ai)r(6 3 - 03) 



,03— ai 



+ 63 - ai - 03) 



x / ^-° 3 - 1 (l - t 2 )^-^-i(l _ s /t 2 ) b i+ b 3-«i-«3-i 2Fl ( 6 3 ai, 61 ai 

\ bx + 6 3 - ax - a 3 



where we used Euler's integral 
1 



l--\dt 2 , (32) 



r b\— 01 — 1 



(1 _ ^fa-oa-i r(6i - oi)r(6 3 - a 3 ) ^ / 61-03,61-01 



2-^1 



(1 + x{t 2 /s - l))6i-«3 r(6i + 63 - 01 - 03) I 61 + 63 - ax - a 3 



1-^ 



and Pfaff 's transformation 

_, / 61 - o 3 , 61 - ai 



2-^1 



61 + 63 - 01 - a 3 

Finally, make the change of variables 

1 - s/t 2 s 
V = — , h 



s , 



1--1 = (t 2 / s yi- bl 2 Fx 



63 - ax, bx - a x 
bx + h-ax- a 3 



1-s 



1-1/(1 -a) 



l-y(l-a) 



(l-y(l-,))S 



in d32j) to get: 

ff((a3); (63)53) 



r(bi -ai)r(&3-a 3 ) _ 6l+62 +63 _ 0l _ a2 _ 03 _ 1 

T(6i + 63 - ai - 03) 1 ; 

1 

y 



X 



bi+b3-ai-o 3 — lM _ y\l>2— a2-l 







(1 - y(l - s)) 6 2-as 



"2-^1 



63 - ai,6i - ax 
bx + 6 3 - 01 - a 3 



1/(1-*) dy. (33) 



This can be further expressed in terms of Appell's function F 3 , defined by \14\ formula 5.7(8)] 



^(ai) "25 /?i,/%; 7; 21,22) = 



(ax) k (a 2 )i(Px)k(f3 2 )i z\ z l 



k,l=0 



(7) 



k+l 



k\ 11 
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Using the double integral representation [TU formula 5.8(3)] 

r(7) 



F 3 (a 1 ,a 2 ;Pi,P2;i;zi,z 2 ) 



r(/3 1 )r(/3 2 )r( 7 -/3 1 -/3 2 ) 



l-u 



(1 - UZl) ai 




(1 - vz 2 ) 



a.2 



dv (34) 



we can obtain by the substitution t = v/(l — u) in the inner integral, an application of the Euler 
integral for 2 F± and changing u — ► y = 1 — it: 



yr-ft-i(i _ y )ft-i / a2i/ g 2 



\ r(7)(i — zi)~ ai 

yz 2 \ dy = j^y^—^jF 3 (a 1 , a 2 ] f3 2 -,Ti zi/(zi -l),z 2 ). 



This formula is a slightly different guise of |26l formula 2.21.1.20]. Combined with (|33p it yields 

n((n,Y (h,Y - r (^l ~ Q l) r (fr3 ~ Q 3 )r(6i + h + fe 3 ~ Ol ~ Q2 ~ Og) q q a +a a -» 1 -6 2 

r(0 2 - a 2 )[r(6i + 63 - ai - a 3 )f 
xtl-^i+fe+fe-ax-as-os-l^^^^ 

(35) 

Finally, we obtain the following representations for 4.F3: 

r(Mr(6 2 )r(& 3 ) 



iF 3 (a, (a 3 ); (63); -z) 



r(ai)r(a 2 )r(a 3 )r(6 2 - a 2 )r(&! + 63 - ai - o 3 ) 



1 1 



s a 2 -lybi+b 3 — ai-a 3 -ln _ s \6i+6 2 +6 3 -ai -02—03-1 H _ y\&2— 02— 1 
X ' ' (l + sz)' T (l-y(l-s)) f> 2- a 3 2i?1 







^3 - aiA - ai 
61 + b 3 - ai - a 3 



2/(1 -s) 



dsdy, 
(36) 



4F 3 (<t, (03)5(63);-*) =-Bi 



s o 2 +03 -6 2 -l/^ _ ^61+62+63-01-02-03-1 
(TTsip 



^3(62 - 03, b 3 - ai;b 2 - a 2 ,h - ar,b 1 + b 2 + b 3 - a x - a 2 - a 3 ; 1 - 1/s, 1 - s)ds, (37) 



where 



£1 



r(&!)r(6 2 )r(6 3 )r(6i + 6 2 + 6 3 - a x - a 2 - a 3 ) 



r( ai )r(a 2 )r(a 3 )[r(6 2 - a 2 )] 2 [r(&! + b 3 - ai - a 3 )] 2 • 

Both representations (|36p and (|37p are believed to be new. They have been verified by termwise 
integration of the series expansions of hypergeometric functions occurring in the integrands and 
comparing coefficients at powers of (—z). According to |26[ formula 7.2.4.74)] the function -F 3 (l — 
1/s, 1 — s) encountered in (|37|) can be expressed as the sum of three 3 F 2 (s). 



5. Open questions and conjectures. The cases q = 1, 2, 3 leave little doubt that the function 
g((a q ); (b q ); s) defined by ([8]) can be expressed in terms of multiple hypergeometric functions for all 
q. The restrictions on parameters in Lemma [1] are rooted in the definition of g as the integral (jSJ). 
These observations motivate 

Open problem 1: How to express the function g((a q ); (b q ); s) defined by (J8]) in terms of 
multivariate hypergeometric functions for q > 3? How to extend the validity of Lemma [1] to a 
wider range of parameter values using the analytic continuation of g((a q ); (b q ); s)? 
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Numerical experiments show that the condition > > is sufficient but in no way necessary 
to the validity of Theorem [TJ Hence, our 

Open problem 2: How to relax the restrictions on parameters in Theorem [T]? 
Numerical tests also indicate clearly the following 

Conjecture 1. Theorem [3] is true for all a > and Yli=i(bi — o-i) > 0. 

The asymptotic formula (|16p and numerical experiments suggest the following 

Conjecture 2. For < a < min(ai,a2, . . . ,a q ), Y^i=xQ } i ~ a i) > an d x > 0: 

q +iF q (a, (a q ); (b q ); -x) < (38) 

Combined with f|28|) Thomae's first relation f|30|) leads to the following curious identity 
i 

y t « (1 _ 4) , + ~ Vl(e _ c , d _ c;d+e _ 6 _ c;1 _ 4)(ii 



T(b)T(c)T(d + e - a - b - c) 
T(a)T(d-a)T(e-a) 



l 



J f-^il-t^-^F^e-^e-fyd + e-b-c-A-^dt. 



o 

Open problem 3: How can one derive the above identity directly from the properties of 2-^1? 
Such derivation would immediately give another proof for the first fundamental relation of Thomae. 
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